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Mean Absolute Value and Standard Deviation of the Phase 
of a Constant Vector Plus a Rayleigh-Distributed 
Vector 

J. Ralph Johler and Lillie C. Walters 

The mean absolute value of the phase and the standard deviation of the phase of a 
constant vector plus a Rayleigh-distributed vector are determined by an evaluation of the 
first and second moment integrals of the probability distribution for various values of average 
relative intensity of the random Rayleigh-distributed component. The results of a quadra- 
ture evaluation of the integrals are tabulated over a wide range of values of average relative 
intensity (& 2 = 0.010 to 1,000). 



1. Introduction to the Problem 

Let E s {t) be a random vector whose amplitude, 
\E s (t)\, is Rayleigh distributed in time and whose 
amplitude, \E s (t)\, and phase angle, Arg E s (t), are 
independently distributed. Rice [l] 1 has introduced 
the concept of a constant vector, E , plus a randomly 
distributed vector, E s (t), to describe electric currents 
or electric fields, E(t)=E J rE s (t), subject to various 
natural fluctuations. Norton has presented tables of 
the cumulative distribution [2] of the amplitude, 
\E(t)\, and its phase [3], 12= Arg E(t), with the relative 
intensity, k 2 =\E r \ 2 /\E \ 2 , as a parameter, where \E T \ 
is the root mean square value of the amplitude \E s (t) | . 
In connection with some recent studies of ionosphere 
roughness, Norton [4] obtained the following proba- 
bility density function 2 of the phase, £>(12), by inte- 
grating over the joint probability distribution given 
by Rice [1] with respect to amplitude, E, the time 
derivative of the amplitude, E' ', and the time deriva- 
tive of the phase, 12': 



where 



'1+JHz exp(2 2 ) [1 + erfO)]} exp(-l/P) (1) 
*=«££ (2 ) 



The phase fluctuations are especially important in 
the performance evaluation of radio navigation and 
phase modulation systems since the behavior of such 
systems is dependent upon the behavior of the argu- 
ment or phase, 12, of the vector, E(t). 

The determination of the mean absolute value, 
|12|, of the phase, 12= arg E(t), and the corresponding 
standard deviation, <m, of a constant vector, plus a 
Rayleigh-distributed vector, E(t), by an evaluation 
of the first and second moment integrals of the proba- 
bility density function, ^(12), for various values of the 
random Rayleigh-distributed component of average 
relative intensity, Jc 2 , is the object of this paper. 

i Figures in brackets indicate the literature references at the end of this paper. 
2 It is of interest to note that H. Bremmer has independently derived the expres- 
sion for p(fi) (private communication). 



2. Theory 

The evaluation of the mean absolute value, |12|, of 
the vector argument or "phase," 12, and the corre- 
sponding variance, era 2 , integrals, 



|G|=2 ri2p(12)dS2, 



and 



co ■ 



2 pl2 2 p(12) d!2, 



(3) 



(4) 



as a function of k 2 is the prime task of this paper. 
The probability density function, p(ti), is symmetric 
about zero, 12=0, and hence the mean, 12, is zero. 
The probability distribution, ^(12), involves the error 
function, 



erf(s)=4= V e~ u2 du } 
V*" Jo 



(5) 



which can be evaluated by the following convergent 
series: 



eTt{Z) VxL l!3 + 2!5 3!7 + * ' J 



(6) 



It is frequently more efficient for large values of z, 
(2>1) to use the asymptotic expansion for real 
values of z, 



V3_ 

(2zV 



l-erf( S )=^exp(- 2 ')[l-^+- 

1-3-5 . 1-3-5-7 "1 { . 

(22 2 ) 3_h (2z 2 Y ' ' "J {7) 

The mean absolute value, |12|, and the variance, 
on, integrals can be represented in the theory of 
Gaussian quadrature [5] as a finite sum, 
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1 r* 

|fl|=-exp (-1/k 2 ) Qpo(Q)dQ 

* JO 



1 M 

=- exp (-1/k 2 ) S W'mOitfbC^+ejr, (8) 

7T W =l 



o-g=i exp {-Ilk 2 ) [ * Sl 2 p (Q)dQ 

K JO 



1 M 

=- exp (-1/k 2 ) S ^Vp (aj+e' M , (9) 

7T TO = 1 

m==l, 2, 3 ... M, 

where e M and eV are error terms which can in general 
be made arbitrarily small by increasing M, 



^ (O)=2xexp (l/k 2 )p(Q), 



(10) 



O w =^[0 6 -OJaw+^[0 6 +Qj, (11) 



and where ft 6 =7r, and 12 a =0, are the upper and lower 
limits of integration, respectively. The weight 
functions can be determined from the limits of 
integration, 



W m =-[n b -n a ]H m 



(12) 



Thus, the integer, M, and the limits of integration, 
O fl , 6 , determine the particular values of the inte- 
grand to be calculated in the quadrature. The 
"universal" constants of the theory of Gaussian 
quadrature, H m , x mj can be determined for any 
given f(x) and various M, 



£ 



M 



f{x)dx=^2 H m f(x m ). 



(13) 



The abscissas, x m , are the roots of the Legendre 
polynomials defined by, 



P (x)-- 


=1 






P:(x) = 


=x 






P>(x) = 


3 2 

=2*- 


i 

2 




Pz(x> 


4^~ 


3 

~2 X 




P<(x) = 


35 . 

T* 1 


15 
4 


* 2 +i 



(14) 



Polynomials of higher degree are determined by use 
of the recursion formula, 

(m+l)P m+1 (x)+mP m ^(x)-(2m+l)xP m (x)=0. 

(15) 



The weight coefficients, H m , can be determined 

2 



from the roots, x m , 



H m = 



{l-xl)[P' m {x n )Y 



(16) 



3I.B.M. 650-407-E.D.P. machine. 

4 The integer to the right of each table entry, if present, indicates the power of 
the factor of 10 by which the number is multiplied, thus positioning the decimal 
point. For example, 5.6513905 -2 means 0.056513905. Note that for & 2 =<». 
<7 Q =-n'/V3 and [G] = »72. 
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Since the available weights, Hm, and abscissas, x my 
are limited [6], (M=48), it is quite possible for very 
precise work to split each integral somewhat arbi- 
trarily but consistent with efficiency, 

f(tt)dQ = f(ti)dQ + f(Q)dQ 

J o a J n a J fid) 

/(Q)dQ+ . . .+ /(Q)dO, (17) 

0(2) jQ(n-D 

where 12 (n) =12 6 , and a specified number of intervals, 
n=l, 2, 3 ... , has been selected with limits of 
integration, fl a , fi (1) , 12 (2) , (3) , . . . Q 6 . Each integral 
is evaluated by the previously described quadrature 
(8, 9) with the abscissas and weights (Af=48). 

3. Computation 

The first and second moment integrals, |12| and 
an 2 , were evaluated between the limits, O a =0, 
12 6 =7T, employing the (M=48) Gaussian abscissas, 
x m , and weights, H m , for values of average relative 
intensity, k 2 , between 0.01 and 1,000 of the random 
Rayleigh distributed component. According to the 
theory of Gaussian quadrature, this integration 
is equivalent to fitting a polynomial of degree 
2m— 1 = 95 at 48 points, to the integrand, which 
points are weighted according to previously described 
rules (12) at the particular values of phase, 12=i2 m (ll). 

The integration was checked at values, A 2 =0.01, 
0.1,1, 10, 100 and 1,000, by splitting the integral into 

two integrals, !2 a =0, fli=w ttb=n, n=2, (17). The 

maximum difference between the value of either |fi| 
or 0-q obtained in this way (17) and that obtained with 
a single integral (9) was ±1 in the eighth significant 
figure. There are indications that the last two 
significant figures could possible be in error in certain 
instances. This is probably due to the electronic 
data processing method 3 rather than the quadra- 
ture. The results of the computation are illustrated 
in figure 1 and are presented in table l 4 . 
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Figure 1. Mean absolute value, |fi|, and standard deviation, <tq, as a function of the average relative intensity, k 2 , of the random 

Rayleigh-distributed component. 



Table 1 



fc2 


|0| 
radians 


radians 


/<2 


|Q| 
radians 


radians 


A;2 


pi 
radians 


radians 


/:2 


!Q| 
radians 


radians 


0.010 
.011 
.012 
.015 


5. 6513905 -2 

5. 9281802 -2 
6. 1929049 -2 

6. 9273763 -2 


7. 0889665 -2 
7. 4368405 -2 

7. 7695630. -2 

8. 6933250 -2 


0.090 
.091 
.092 
.095 


1. 7204794 -1 
1. 7303579 -1 
1. 7401916 -1 
1. 7694134 -1 


2. 1758636 -1 
2. 1886360 -1 
2. 2013540 -1 
2. 2391945 -1 


0.80 
.81 

.82 
.85 


5. 7651767 -1 
5. 8018927 -1 
5. 8382359 -1 
5. 9450880 -1 


7.8540S00 1 
7.9022010 1 
7.9497130 -1 
8. 0887390 -1 


7.0 
7. 1 

7.2 
7.5 


1. 1573986 
1. 1601429 
1. 1628333 
1. 1705980 


1. 4351373 
1.4378157 
1. 4404390 
1. 4479959 


.020 
.021 
.022 
.025 


8. 0059192 -2 
8. 2050463 -2 
8. 3996053 -2 
8. 9586354 -2 


1. 0051253 -1 
1.0302171 -1 
1. 0547390 -1 
1. 1252438 -1 


.10 
.11 
.12 

.15 


1. 8172472 -1 
1.9099601 -1 

1. 9992282 -1 

2. 2508374 -1 


2. 3012770 -1 
2. 4221559 -1 
2. 5393404 -1 
2. 8749828 -1 


.90 
.91 

.92 
.95 


6. 1162274 -1 
6. 1494579 -1 
6. 1823668 -1 
6. 2792124 -1 


8. 3093795 -1 
8. 3519350 -1 
8. 3939890 -1 
8. 5172250 -1 


8.0 
8.1 
8.2 
8.5 


1. 1826098 
1. 1848849 
1. 1871207 
1. 1936018 


1. 4596450 
1. 4618459 
1. 4640069 
1. 4702620 


.030 
.031 
.032 
.035 


9. 8223130 -2 
9. 9864166 -2 
1. 0148000 -1 
1. 0618661 -1 


1. 2342886 -1 
1. 2550275 -1 
1. 2754525 -1 
1. 3349834 -1 


.20 
.21 
.22 
.25 


2. 6329686 -1 
2. 7053872 -1 
2. 7767192 -1 
2. 9848460 -1 


3. 4032320 -1 
3. 5059207 -1 
3. 6077869 -1 
3. 9085106 -1 


1.0 
1.1 
1.2 
1.5 


6. 4346186 -1 
6. 7247370 -1 

6. 9902834 -1 

7. 6682649 -1 


8. 7133750 -1 

9. 0744135 -1 
9. 3992080 -1 
1. 0205470 


9.0 
9.1 
9.2 
9.5 


1. 2037103 
1. 2056364 
1. 2075323 
1. 2130479 


1. 4799901 
1. 4818398 
1. 4836595 
1. 4889462 


.040 
.041 
.042 
.045 


1. 1361951 -1 
1. 1505149 -1 
1. 1046690 -1 
1. 2(161991 -1 


1. 4291103 -1 
1. 4472622 -1 
1. 4652096 -1 
1. 5179011 -1 


.30 
.31 
.32 
.35 


3. 3147043 -1 
3. 3783684 -1 
3. 4413045 -1 
3. 6258665 -1 


4. 3923906 -1 
4. 4863226 -1 
4. 5792448 -1 
4. 8517326 -1 


2.0 
2.1 
2.2 
2.5 


8. 5195625 -1 
8. 6599590 -1 
8. 7925345 -1 
9. 1498336 -1 


1. 1175455 
1. 1331300 
1. 1477455 
1. 1866651 


10 
11 
12 
15 


1. 2217077 
1. 2372944 
1. 2509646 
1. 2836835 


1. 4972268 
1. 5120705 
1. 5250259 
1. 5558026 


.050 
.051 
.052 
.055 


1. 2726013 -1 
1 . 2854991 -1 
1 . 2982785 -1 
1. 3359414 -1 


1. 6022582 -1 
1. 6186600 -1 
1. 6349168 -1 
1. 6828606 -1 


.40 
.41 

.42 
.45 


3. 9199002 -1 

3. 9767434 -1 

4. 0329502 -1 
4. 1978140 -1 


5. 2841280 -1 
5. 3672815 -1 
5. 4493255 -1 
5. 6888385 -1 


3.0 
3.1 
3.2 
3.5 


9. 6399295 -1 
9. 7254727 -1 
9. 8075157 -1 
1. 0034854 


1. 2389914 
1. 2480052 
1. 2566181 
1. 2803230 


20 
21 
22 
25 


1. 3212337 
1. 3271198 
1. 3326079 
1. 3470845 


1. 5907359 
1. 5961754 
1. 6012386 
1. 6145548 


.060 
.061 
.062 
.065 


1. 3966386 -1 
1. 4084930 -1 
1.4202581 -1 
1. 4550327 -1 


1. 7602361 -1 
1. 7753642 -1 
1. 7903826 -1 
1. 8348099 -1 


.50 
.51 
.52 

.55 


4. 4604899 -1 
4. 5112729 -1 
4. 5614885 -1 
4. 7088003 -1 


6. 0664350 -1 
6. 1388105 -1 
6. 2101680 -1 
6. 4182710 -1 


4.0 
4.1 
4.2 
4.5 


1. 0361770 
1. 0420645 
1. 0477624 
1. 0638134 


1. 3140093 
1. 3200272 
1. 3258374 
1. 3421326 


30 
31 
32 
35 


1. 3662719 
1. 3695492 
1. 3726740 
1. 3812409 


1. 6321174 
1. 6351074 
1. 6379555 
1. 6457515 


.070 
.071 

.072 
.075 


1. 5113924 -1 
1. 5224405 -1 
1.5334169 -1 

1. 5659419 -1 


1. 9069203 -1 
1. 9210734 -1 
1. 9351423 -1 
1. 9768590 -1 


.60 
.61 
.62 
.65 


4. 9436251 -1 

4. 9890499 -1 

5. 0339766 -1 
5. 1658413 -1 


6. 7460915 -1 
6. 8089385 -1 
6. 8709150 -1 
7. 0517715 -1 


5.0 
5.1 

5.2 
5.5 


1. 0875730 
1. 0919370 
1. 0961852 
1. 1082832 


1. 3660625 
1. 3704336 
1. 3746814 
1. 3867404 


40 
41 

42 
45 


1. 3933403 
1. 3954936 
1. 3975702 
1. 4033835 


1. 6567297 
1. 6586796 
1. 6605589 
1. 6658143 


.080 
.081 
.082 
.085 


1. 6188698 -1 
1.0292752 -1 
1. 6396238 -1 
1. 6703356 -1 


2. 0448601 -1 
2. 0582473 -1 
2. 0715680 -1 
2. 1111367 -1 


.70 
.71 
.72 

.75 


5. 3763008 -1 
5. 4170535 -1 
5. 4573767 -1 
5. 5758218 -1 


7. 3371715 -1 
7. 3919745 -1 
7. 4460535 -1 
7. 6040620 -1 


6.0 
6.1 
6.2 
8.5 


1. 1265439 
1. 1299436 
1. 1332664 
1. 1428021 


1. 4048374 
1. 4081929 
1. 4114684 
1. 4208462 


50 
51 
52 
55 


1. 4118941 
1. 4134456 
1. 4149523 
1. 4192256 


1. 6734927 
1. 6748905 
1. 6762476 
1. 6800932 
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Table 1. — Continued. 



fc2 


|0| 
radians 


radians 


W 


radians 


radians 


60 
61 
62 
65 


1. 4256277 
1. 4268134 
1. 4279704 
1. 4312811 


1. 6858461 
1. 6869104 
1. 6879488 
1. 6909181 


400 
410 
420 
450 


1. 5143646 
1. 5150558 
1. 5157220 
1. 5175857 


1. 7646036 
1. 7652101 
1. 7657947 
1. 7674297 


70 
71 
72 
75 


1. 4363217 
1. 4372656 
1. 4381902 
1. 4408523 


1. 6954337 
1. 6962787 
1. 6971060 
1. 6994873 


500 
510 
520 
550 


1. 5203113 
1. 5208079 
1. 5212898 
1. 5226565 


1. 7698191 
1. 7702541 
1. 7706766 
1. 7718738 


80 
81 

82 
85 


1. 4449535 
1. 4457283 
1. 4464888 
1. 4486897 


1. 7031527 
1. 7038446 
1. 7045237 
1. 7064882 


600 
610 
620 
650 


1. 5247021 
1. 5250810 
1. 5254503 
1. 5265072 


1. 7736652 
1. 7739967 
1. 7743202 
1. 7752451 


90 
91 
92 
95 


1.4521112 
1. 4527618 
1. 4534017 
1. 4552605 


1. 7095402 
1. 7101201 
1. 7106904 
1. 7123466 


700 
710 
720 

750 


1. 5281154 
1. 5284164 
1. 5287113 
1. 5295601 


1. 7766522 
1. 7769155 
1. 7771734 
1. 7779157 


100 
110 
120 
150 


1.4581716 
1. 4633887 
1. 4679420 
1. 4787613 


1. 7149389 
1. 7195798 
1. 7236250 
1. 7332182 


800 
810 
820 
850 


1. 5308671 
1. 5311138 
1. 5313563 
1. 5320572 


1. 7790585 
1. 7792741 
1. 7794859 
1. 7800987 


200 
210 
220 
250 


1. 4910552 
1. 4929717 
1. 4947562 
1. 4994527 


1. 7440872 
1. 7457784 
1. 7473525 
1. 7514919 


900 
910 
920 
950 


1. 5331466 
1. 5333538 
1. 5335573 
1. 5341489 


1. 7810506 
1. 7812315 
1. 7814094 
1. 7819263 


300 
310 
320 
350 


1. 5056546 
1. 5067115 
1. 5077186 
1. 5104767 


1. 7569509 
1. 7578805 
1. 7587659 
1. 7611897 


1000 


1. 5350752 


1. 7827352 



4. Discussion and Conclusions 

The problem formulated in this paper was the 
evaluation of two integrals — the first and second 
moments of the probability distribution of the phase 
of a constant vector plus a Kayleigh-distributed 
vector. These integrals were evaluated by an 
application of the theory of Gaussian quadrature and 
the results indicate that the method provides an 
efficient and precise solution to the problem. 



It is quite possible to estimate the entire proba- 
bility distribution from physical measurements of 
the standard deviation or the mean absolute value 
of the phase with the aid of the table or figure. 
This work therefore provides an analysis tool for the 
study of phase fluctuations in radio navigation and 
phase modulation systems. 



This work was suggested by K. A. Norton and 
performed under National Bureau of Standards 
Project 8330-11-8332. The careful electronic com- 
puter work was developed by L. C. Walters. 
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